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The dynamics of quantum degenerate two-dimensional dipolar excitons confined in electrostatic 
traps is analyzed and compared to recent experiments. The model results stress the importance 
of artificial trapping for achieving and sustaining a quantum degenerate exciton fluid in such sys- 
tems and suggest that a long-lived, spatially uniform, and highly degenerate exciton system was 
experimentally produced in those electrostatic traps. 



Excitons are excellent candidates for producing Bose- 
Einstein (BE) quantum statistical effects, such as con- 
densation and superfluidity, on a semiconductor platform 
P, . Observation of such effects can open up new and 
exciting opportunities for both fundamental and applied 
research. A crucial step in this direction is the ability 
to obtain long-lived quantum degenerate excitons, which 
practically means sustaining a high density exciton cloud 
for times long enough for the excitons to cool down to 
the lattice temperature. One of the more promising two- 
dimensional (2D) exciton systems for this purpose is the 
dipolar (or "indirect") exciton system in double quan- 
tum well (DQW) structures. In this system, the con- 
stituent electrons and holes reside in difference quantum 
wells, giving the excitons very long recombination life- 
times (in the microsecond range) which should be much 
longer than their thermalization time with the host lat- 
tice. In addition to being able to cool the excitons, to 
observe BE effects, one must also be able to obtain high 
densities. Unfortunately, the combined task of cooling an 
exciton fluid while maintaining high density has turned 
out to be quite difficult, as recent studies have revealed 
HQ: the inherent strong dipolar repulsion between pairs 
of excitons, while desirable for eliminating the forma- 
tion of eh droplets and other complexes, also leads to 
a rapid spatial expansion of the initially dense exciton 
cloud. As a result, the exciton density quickly drops to 
below the critical value for quantum degeneracy. In addi- 
tion, the strong exciton expansion prevents the formation 
of a static fluid and the dipolar repulsion energy is con- 
verted into exciton kinetic energy during the expansion, 
which can result in a heating source which slows the cool- 
ing process. This fast expansion of the exciton cloud thus 
poses a serious technical problem to producing BE effects 
in exciton systems. 

One possible solution is to continuously supply cold 
dipolar excitons, using long distance in-plane charge 
transport ("exciton rings") as a cooling mechanism jj, 
E| Q- Recently we have been trying to find a solution 
by implementing spatially confining potential traps for 
the dipolar excitons in the plane of the quantum wells, 
an approach that is somewhat analogous to that of atom 
trapping and cooling. Two different trapping methods 
have been attempted exp er iment all y : strain-induced || 
and electrostatic |slllfj,lll| traps (Xtraps). The Xtraps, 



the design of which is given in Ref. Q seem to be par- 
ticularly promising as they are easily tunable in size and 
depth and have sharp boundaries, all of which are impor- 
tant for achieving high exciton densities, as was recently 
demonstrated experimentally [lfj. 

In light of these recent experiments, it is important 
to analyze and understand the exciton dynamics in such 
Xtraps and their critical role in achieving a long life- 
time and spatially uniform degenerate fluid of excitons. 
The calculations in this paper are based on the models 
developed by Ivanov et. al. and extended to in- 

clude the Xtrap potential and expansion induced exciton 
heating. They explain our previous experimental obser- 
vations and reveal the detailed dynamics of the dipolar 
exciton fluid confined in a trap. In particular, these cal- 
culations show the advantages of this trapping method, 
and suggest that the exciton fluid confined in the Xtraps 
is highly degenerate and spatially uniform, over long time 
periods of the order of the exciton lifetime. 

The inset of Fig. ^ shows an electrostatic trap design 
0, where dipolar excitons (with a dipole moment given 

by dx = —ezQZ, zq being the effective dipole length) are 
trapped under a local circular electrostatic gate with ra- 
dius R. When an electrical bias is applied to the gate, 
the (negative) dipole-field interaction energy profile et at 
the DQW position, plotted in Fig. (solid line), effec- 
tively confines the excitons to the region under the cir- 
cular gate with the trapping force given by Ft = — Vet 
|9|. For calculation convenience, we approximate (Fig.^ 
dashed line) the xtrap potential with an analytic func- 

tanh 2 ^ r 7 R ^ ) . Here, e° is the energy 
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tion: e t (r) = -f ( 1 



at the trap center, which is proportional to the applied 
voltage, and 5 is the effective " thickness" of the bound- 
ary region, proportional to the distance between the the 
top gate and bottom electrode. 

Two scattering mechanism for dipolar excitons are in- 
cluded in this model. The first one is the fast, density 
dependent exciton-exciton scattering, which is dominant 
at high exciton densities. Such a scattering process, with 
a characteristic time much smaller than any other time 
scale of the dynamics, yields internal exciton equilibrium 
and therefore a well defined exciton temperature, Tx- It 
is also responsible for a density dependent diffusion co- 
efficient, Dxx, calculated by Ivanov et al. in Ref. fl2j |. 
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FIG. 1: The calculated Xtrap potential radial profile (solid 
line) and its analytic approximation (dashed line). The inset 
shows an schematic diagram of an Xtrap design geometry. 



The other scattering mechanism is due to exciton-QW 
disorder interaction, which is density independent but 
yields a diffusion coefficient, Ddo, which depends on the 
QW width to the sixth power (see Refs. 0, 0) We 
also assume that there is a single characteristic time (ji ) 
for the thermalization of the exciton gas with the lat- 
tice (with a lattice temperature Ti). The dynamics of a 
dipolar exciton fluid in an Xtrap can then be described 
by two coupled equations for the exciton fluid density 
and temperature. The first nonlinear diffusion equation 
describes the time evolution of the exciton fluid: 
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Ix(r,t)=0. (1) 



Here, nx = nx(r,t) is the exciton density profile. 3d, 
3d, and Jt are the three exciton currents which are driven 
by diffusion, dipolar repulsion, and the trapping poten- 
tial, respectively. Ix(r,t) = Ix(r)S(t) is the exciton 
source for a short optical pulse excitation and Tx is the 
lifetime of the excitons, which is assumed to be density 
independent. These exciton currents are related to the 
various forces through the exciton mobility J = nxnF. 
The exciton mobility, a, is related to the effective ex- 
citon diffusion coefficient D — DxxDd /(Dxx + Ddo), 
through Einstein relations u = (D/kTo)(e T °/ Tx — 1), with 
To = (2irh 2 nx) / (kgmx) being the degeneracy tempera- 
ture. Here, mx = 0.2m e is the exciton mass and g = 4 
is the exciton spin degeneracy. The dipole-repulsion 
force, Fd = —aVnx is due to the mean field dipole- 
dipole interaction energy 01 : Edd = iiredxnx I ' e = oinx, 
where e is the background dielectric constant. The dif- 



fusive force is 
C = kT x ln (l - 
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jVnx, where 



is the chemical potential in the 



non-interacting limit. The trapping force for the dipo- 
lar excitons can be derived from the above trapping po- 
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pressions for the exciton currents in Eq. are there- 



fore: 
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The second equation describes the time evolution of 
the exciton fluid temperature, Tx- 
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where Edd, the total potential energy due to the dipolc- 
dipole interactions is given by: Edd = f an x (r,t)d 2 r 
and Nx(t) = J nx(?,t)d 2 r is the total number of dipolar 
excitons. The first term on the RHS of Eq. [3 represents 
the heating of the exciton gas due to the driven expansion 
and the last term represents the thermalization of the 
excitons with the lattice. 

Eqs.Q,© are coupled and can be numerically evalu- 
ated. For the following calculations, we use parameters 
that are typical for our experiments with GaAs double 
QW structures [T^- The size of the trap is taken as 
R = 25 um, consistent with the experiments performed 
in ^3] . The depth of the trap is fixed at e° = 30 meV. 
For a short (~ Ins) pulse excitation with a gaussian pro- 
file, we assume that dipolar excitons are subsequently 
created at the center of the trap with an initial gaussian 

2/2 

profile: nx{r, t = 0) = n§e~ r > w ° where typically for our 
experiments, wq — 15 um and no = 1.5 x 10 11 cm -2 . 
The initial exciton temperature is taken to be the energy 
difference between the optical excitation energy and the 
dipolar exciton energy, Tx(t — 0) ~ 50^, which is much 
hotter than the lattice temperature T; = 1AK, as we 
assume that the optically excited carriers rapidly trans- 
form into hot dipolar excitons. The exciton lifetime, tx 
is position dependent. Inside the trap, the excitons are 
dipolar and their lifetime is in microsecond range. For 
the particular trap depth 30 meV, we use a typical value 
of tx = t~x id = 1.5//S. Outside the trap, the excitons are 
direct and their lifetime is taken as tx = Tx n = 0.1ns. 
In the trap boundary region, the dipolar excitons are 
subject to ionization due to the in-plane fringing fields 
that depends on the geometry of the trap design (see 
Ref. |9( for details). This ionization process yields an 
effective boundary lifetime for the excitons (r trap ). We 
have shown analytically 0] that this lifetime, and thus 
the right Xtrap design, is crucial for the effectiveness of 
the Xtrap, as was further confirmed by our experiments 
|To| and the following calculations. 

Fig-EK plots the radial density profile of the excitons in 
a high quality Xtrap with negligible boundary ionization, 
in which case the boundary lifetime is only limited by the 
intrinsic lifetime of the excitons (r trap = Tx)- The val- 
ues of the diffusion coefficients are chosen by fitting the 
expansion dynamics (the cloud FWHM as a function of 
time in particular) of a free expanding excitons to ex- 
perimental results (see Ref. 0). In a good agreement 
with the experimental results of Ref. [10J, within tens 
of nanoseconds the profile flattens, as the driven expan- 
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FIG. 2: Radial profile of the exciton density distribution for 
(a) a high quality Xtrap (jtrap = 1.5/is) at various times, 
and (b) High quality Xtrap (rtrap = 1.5/xs, solid line), low 
quality Xtrap (jtrap = 0.1ns, dashed line), and free expanding 
excitons (dotted line) at t = 45ns. The bar at the bottom of 
each figure mark the radial extent of the Xtrap. 
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FIG. 3: The exciton ground state occupancy, Ne=o, as a 
function of time for a high quality Xtrap (jtrap = 1.5us, solid 
line), a low quality trap (rtrap = Ins, dotted line) and free 
expanding excitons (dashed line). The inset shows the cor- 
responding exciton temperature as a function of time for the 
high quality Xtrap (solid line) and free expending excitons 
(dotted-dashed line). 



sion by the dipole repulsion pushes the excitons to the 
reflecting trap boundary. These two competing forces re- 
sult in a stable flat exciton density profile (the density 
profile with the lowest total energy) that decays with the 
characteristic radiative lifetime of the excitons, tx- As 
we will show below, as the profile flattens out, the exci- 
ton quickly cools to the lattice temperature and becomes 
highly degenerate. This degeneracy is maintained on the 
time scale of the exciton radiative lifetime. 

This result is strikingly different from either a free ex- 
panding exciton cloud or excitons in a low quality Xtrap 
with a considerable boundary ionization, as depicted in 
Fig. 12b. In the case of free expanding excitons, the den- 
sity rapidly decreases, as the driven expansion contin- 
ues to spread the exciton cloud to large radii. For a 
low quality trap (Tt rap = 0.1ns), excitons reaching the 
trap boundaries are effectively eliminated by ionization, 
rapidly depleting the expanding exciton cloud. Note that 
in this case, the exciton profile remains curved and does 
not flattens like the exciton pool in the high quality trap. 
This is also in agreement with the experimental profile 
found for a low quality Xtrap in Ref. [lJJ. 

To see the dramatic effect of trapping on high mo- 
bility excitons (Ddo — 10cm 2 /s), we show the calcu- 
lated ground state occupation number of the 2D excitons, 
Ne=o = (exp(Xb/7x) — 1), as a function of time after 
photoexcitation, for excitons in a high and low quality 
Xtrap as well as for free excitons. Initially, the ground 
state occupation is much smaller than unity, as the exci- 
tons are hot and essentially classical. The exciton cloud 
then start to expand and cool. As the transition to a 
distinctive BE statistics arises when T /Tx ~ 1, there 
is a competition between cooling, that tend to increase 
Ne=o and expansion that tend to decrease it due to the 
drop of nx- In a high quality Xtrap, expansion is lim- 



ited to the Xtrap boundaries, and thus reaches steady 
state after ~ 10 nanoseconds. Cooling is then efficient 
(as is seen in the inset of Fig. [3J and the excitons reach 
a stable high degenerate state with Ne=o ^> 1, that de- 
cays slowly on the time scale of tx ■ On the contrary, if 
the excitons are freely expanding (no trap) , their density 
continues to drop as they expand. While they continue to 
cool by interaction with the lattice, the expansion itself 
heats up the cloud, as more internal potential energy is 
lost and converted into heat by the fast exciton-exciton 
interaction. This leads to a slower net cooling rate for the 
free expanding excitons compared with the trapped ones, 
and consequently the excitons reach only marginal quan- 
tum degeneracy before their density drops too low due to 
the fast driven expansion. Finally, in a low quality trap, 
the drop in exciton density due to boundary ionization 
dominates over the cooling, completely preventing the 
excitons from reaching degeneracy. 

While those calculations strongly suggest that the exci- 
ton pool that was created in Ref. [10( is highly degenerate 
for hundreds of nanoseconds, it is based on the assump- 
tion that the cooling rate is indeed much shorter than 
the exciton lifetime. While this is a very reasonable as- 
sumption (see Ref. [I2|)j a clear experimental evidence 
of this degeneracy that is independent on model assump- 
tions is still essential, and experiments looking for such 
conclusive evidence are ongoing. 

In summary, the dynamics of two-dimensional dipolar 
excitons confined in electrostatic traps is modeled and 
agrees well with recent experiments. The model shows 
that artificial trapping is crucial for achieving and sus- 
taining a quantum degenerate exciton fluid is such sys- 
tems and suggest that a spatially uniform, highly degen- 



4 



erate exciton pool was experimentally produced in those 
electrostatic traps. 
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